Dipolar ground state of planar spins on triangular lattices 
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An infinite triangular lattice of classical dipolar spins is usually considered to have a ferromagnetic 
ground state. We examine the validity of this statement for finite lattices and in the limit of large 
lattices. We find that the ground state of rectangular arrays is strongly dependent on size and 
aspect ratio. Three results emerge that are significant for understanding the ground state properties: 
i) formation of domain walls is energetically favored for aspect ratios below a critical value; ii) a 
vortex state is energetically favored in the thermodynamic limit of an infinite number of spins, but 
nevertheless such a configuration may not be observed even in very large lattices if the aspect ratio 
is large; iii) finite range (R) approximations to actual dipole sums may give spurious results and the 
limit R — > oo depends on the way it is taken. For the usual, isotropic limit, the ferromagnetic state 
is linearly unstable and the domain wall energy is negative for any finite range cutoff. 

PACS numbers: 75.10.Hk, 75.75.+a, 75.60.Ch, 77.80.-e 
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Electric or magnetic fields generated by point dipole 
sources have significant effects on polarizable materials. 
The effects can dominate in thin film and nanostructured 
elements. In one and two dimensions, long range order 
of electric or magnetic polarizations in dipole lattices is 
possible! and it is particularly interesting in view of the 
very strong impact that lattice geometry can have on 
the nature of the long range order £ Interest in these ef- 
fects has been renewed by the ability to create arrays 
out of polarizable material a 3 ^ 5 with element sizes that 
can behave to reasonable approximation as point dipole 
sources. 

Whereas long range ordering in infinite two dimen- 
sional lattices has been examined in some detail for differ- 
ent rhombic lattices*^ there is to date very little known 
about how exactly long range order is affected by lat- 
tice size, shape and dimension in the transition region 
between small and infinite lattices. We find that this is 
a non-trivial point for two dimensional lattices, and note 
that many real systems of technological interest fall in 
the very large, yet non-infinite, category. Our results im- 
pact not only the nature of ordering in large, rectangular- 
lattices, but also some inadequacies and failures in the 
use of an interaction cutoff for the dipolar interaction 
(see also Refs.|H§). 

In this manuscript we study a two-dimensional triangu- 
lar lattice of classical planar spins with long range dipolar 
coupling between spins. The nearest neighbour distance 
between spins is taken as unit length. The total dipolar 
energy is 
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, (Si ■ fij)(Sj ■ Tij) 



(1) 



where (ij) is a couple of spins and is their distance. 
The energy can also be written in terms of an effective 
field Hi acting at a lattice site i as E = — § J2i ' Hi- 
Each Si may represent either a single spin or the effec- 



tive spin of a single domain magnetic particle. It has 
been showniS that near-field corrections due to the finite 
size of magnetic particles can be accounted for with cor- 
rection terms that do not alter the general qualitative 
features of the planar array interactions. The prefactor 
to = \9' ''^b^ 2 1S positive and absorbs the modulus of 
the spin (so that \Si\ — 1). The prefactor can be set to 
51 = 1 , if there are no other terms in the energy. 

The dipolar interaction between two spins S± and S2 
located at sites r\ and ?% is minimized by aligning them 
ferromagnetically along the line fx% = f± — r*2- A sin- 
gle chain of spins will therefore order ferromagnetically 
along the chain. For parallel chains of spins, the field 
H(d) generated by an infinite chain at a point some dis- 
tance d away from the chain axis plays a special role be- 
cause H (d) = in the continuum approximation^ For 
a discrete lattice, H(d) does not vanish, but it decays 
exponentially 1 !^ with the distance, with a sign depend- 
ing on the precise discrete structure. 

As a consequence of thatji 2 - geometry becomes espe- 
cially important in two dimensions: the sign of H(d) de- 
pends on the lattice structure and is negative for a square 
lattice but positive for a triangular latticed Spin chains 
are ferromagnetically coupled in a triangular lattice and 
antiferromagnetically coupled in a square lattice. 

This picture breaks down when the chains are finite in 
length. Consider a rectangular sample of sides L x and 
L y with aspect ratio r = L x /L y > 1. We shall show, 
using direct numerical evaluation of Eq. ft), that the 
energy of a domain wall parallel to the x axis scales with 
L x and becomes negative for r < r c where the critical 
aspect ratio is r c w 4.8. This implies that domain wall 
formation is favored in rectangles with relatively small 
aspect ratios. Most significantly, increasing the size of 
the sample while keeping r fixed, does not increase the 
number of domain walls. Whereas one might anticipate 
that array shape can drive domain formation as it does in 
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exchange coupled ferromagnets, the number of domains 
in a dipole ferromagnet is size independent. This feature 
is not found in exchange coupled ferromagnets. The dif- 
ference is because the non-local dipole interaction in a 
dipole array produces both the magnetic anisotropy and 
the magnetic coupling required for formation of a domain 
boundary wall. In the exchange coupled ferromagnet, ex- 
change interactions are usually negligible beyond nearest 
or next-nearest neighbors. 

In the following we shall discuss these shape and size 
effects, and also show that vortex formation is energeti- 
cally favored only for very large lattices when r is large. 
Appearance of vortices is a well known fact, both in 
purely dipolar systemsA&i^ and in dipole plus exchange 
systems^ but we are able to evaluate the vortex energy 
as a function of size and shape of the sample. Effects 
of local quadratic and quartic magnetic anisotropics are 
also discussed. 

We begin by reviewing two results which are particu- 
larly enlightening as to why ferromagnetic lines and vor- 
tices are useful concepts for understanding ordering in 
two-dimensional dipolar lattices. Let us consider a closed 
curve r. If the spin is treated as a continuous field, we 
have 5(r) at each point f of the curve. The first result is 
that, if S is tangent to T, the dipolar field H(P) gener- 
ated in any point P not belonging to the curve vanishes. 16 
This result shows that concentric circles in a vortex or 
circles belonging to neighbouring vortices do not inter- 
act for continuous dipole sources. If the discrete lattice 
structure is considered, this is no more exact, but the in- 
teractions are weak, specially for non concentric circles, 
as it occurs in neighbouring vortices. 

The second result is that the field produced by all spins 
of a closed curve T in a point P belonging to the curve 
is tangent to the curve itself>i£ Therefore the vortex is 
the spin configuration that minimizes the dipolar energy 
of a curve of spins, as the ferromagnetic configuration 
minimizes the dipolar energy of a straight line of spins. 

Two conclusions follow. First, ferromagnetic lines or 
vortices are energetically favorable because in such con- 
figurations each spin is aligned along the dipolar field. 
Second, interaction between (infinite) straight lines or 
vortex lines is weak (in a continuum approximation it 
is exactly zero) and exists only in discrete lattices. The 
type of alignment therefore depends not only on size and 
shape of an array, but also on the local symmetry defining 
the geometry of the array. 

Now let us consider a rectangular sample of sides 
L x ,L y with L x > L y . Because of shape effects^ the 
ferromagnetic state is directed along the x axis. We cre- 
ate a domain wall along the x axis and form two equally 
large domains above and below it, with the magnetiza- 
tion pointing in ±x respectively. In the inset to Fig. ^ we 
compare the energy of the two-domain state (E2) with 
the single-domain (FM) state. Shown is the quantity 
Edw = (E2 — Efm)/L x which is a function of the as- 
pect ratio r only. The formation of a domain wall defin- 
ing the two-domain state is energetically favored when 
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FIG. 1: Inset: difference between the energy of the two- 
domain state, E2, and the energy of the single-domain state, 
Efm, per unit length L x of a rectangular sample of sides 
L x , L y , as a, function of the aspect ratio r — L x /L y . The 
plotted quantity is also called domain wall energy, Edw Main: 
energy of a vortex (Ev) as a function of the horizontal size 
L x . Again, we scale with respect to Efm and per unit length. 
Straight lines are fits to the equation (Ey — Efm)/L x = 
a(r) - b(r) \nL x . 



r < r c fts 4.8. For larger r, the formation of the domain 
wall is energetically unfavorable. A distinguishing fea- 
ture of the dipole lattice is that the energy of the domain 
wall per unit length depends on the aspect ratio only, not 
on the size of the sample. One does not see a proliferation 
of domains with increasing array size. 

We now examine the energy of a vortex state in a rect- 
angular sample. For the sake of simplicity, we consider 
an elliptic vortex, as defined in note fl9l This config- 
uration is not expected to be the 'exact' ground state, 
because of corner effects. A rigorous treatment should 
allow for spin relaxation, which has been done in a few 
cases (see below). However, it is not feasible to do that 
for a large number of big samples of different aspect ra- 
tios. We are confident that our main results, Eq. © and 
Figs. I (main) and 2, are robust, because corner effects 
should be negligible in the limit of big sizes. 

In Fig. n (main) we show the energy Ey of a vortex 
as a function of the size L x , for different aspect ratios 
r. Unlike the domain wall energy, which depends only 
upon r, the quantity (Ey — Efm)/L x depends both on 
the shape and the size of the sample. We find that the 
vortex energy obeys the functional form 



{Ey - E FM )/L X = a{r) - b(r) lnL a 



(2) 



with a(r) and b(r) plotted in Fig. El The first term, a(r), 
does not depend on L x and displays an r-dependence 
which strongly resembles Ed w (r) (see the inset of Fig. 
An average of the magnetization of the vortex in the up- 
per (y > 0) and lower (y < 0) regions approximates a 
two-domain state, and a(r) describes the interaction be- 
tween these upper and lower parts of the vortex. 

The second term, — b(r)hiL x , is always negative and 
drives vortex formation. Consider a two-domain ferro- 
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FIG. 2: Dependence of the parameters a(r) (open circles) and 
b(r) (open diamonds) on the aspect ratio r of a rectangular 
sample. Their meaning is illustrated in Fig. Q (main). 



magnetic alignment as an initial condition on the array. 
The energy gain for deformation into a vortex increases as 
L x is increased, but strongly decreases as r is increased. 
Therefore, provided that the size L x is sufficiently large, 
the vortex state appears to have an energy lower than 
parallel aligned states for any aspect ratio. The size nec- 
essary for this depends on 

Our main study is based on the direct numerical eval- 
uations of dipolar energy for a number of fixed ini- 
tial spin arrangements. This method allows to deter- 
mine the exact energies for specified configurations in 
very large lattices. In order to take into account the 
possibility of spin relaxation, in a few cases configura- 
tions that produce stable configurations have also been 
found using numerical methods that adjust local spin ori- 
entations in such a way as to minimize the total energy 
iteratively. The results are generally consistent with the 
main study: Vortex-like configurations appear for large 
enough arrays (large L x ), and can evolve into some type 
of domain structure that reduces magnetostatic energy 
near the array edges. For large values of the aspect ra- 
tio vortices do not form, but domains form in proximity 
of the smallest sides because a FM state would have the 
local magnetization perpendicular to them, which is en- 
ergetically unfavourable. These edge effects do not seem 
to penetrate the whole sample for large enough arrays. 
The actual configuration is determined by instabilities of 
parallel rectangular domains with respect to the forma- 
tion of flux-closure domains with triangular shape. More 
general and precise conclusions would require a system- 
atic simulation work^i which is beyond the scope of this 
Rapid Note. 

We now turn to approximations involving the range R 
of the dipolar interaction. Truncation of dipole sums 9 
by imposing finite range cutoff's for dipolar interactions 
is often used in numerical simulations in order to reduce 
computation times. This may be a severe approximation 
for dipole lattices regardless of lattice size. 

The stability of a ground state configuration can be 
examined by studying the frequencies of spin wave exci- 



tations. It is widely accepted that, for an infinite range 
interaction, the ferromagnetic state on a triangular lat- 
tice is locally stable. In consequence, the linear spin wave 
frequency spectrum is real over the entire Brillouin zone, 
and co(q) « ^fq at small q. The non-analyticity of Lo(q) 
in q = is due to the long-range character of the inter- 
action. 

An estimate of the general behavior of spin wave fre- 
quencies in parallel (FM) spin configuration under the 
finite range approximation can be made as follows. If the 
spin configuration is assumed to consist of infinite par- 
allel lines it is easy to calculate the domain wall energy 
because interaction between lines decays exponentially, 
and we get E^ = 0.2356. This positive value means 
that creating a domain wall in an infinite system with 
infinite range interaction has a finite energy cost. Now 
let the range R be finite. 

The spin wav e frequency has the following general 
formi 2 ^ w(q) = y/(A — B)(A + B), where 



A - B = 



1 



[2 + cos(g • n)} 



A + B = D yy (q) ~ D yy (0) 



(3) 
(4) 



with Dyy(q) = J2i y- ~ 3fjrJ cos(g • n). 
The quantity (A — B) is manifestly and strictly positive, 
while (A + B) may take both signs, as we are going to ar- 
gue. For infinite range interactions^ D yy (q) — D yy (0) w 
cq, with a positive prefactor c. For a finite range R and 
small gr, we can expand the cosine in D yy and replace the 
discrete summation on the angle 9 [fi = (r;,#i)] with an 
integral, getting 



2 

A + B^lY — (6sin 2 0-l) 



(5) 



where </> is the angle between q and the x axis. 

The above summation is limited by the range R of the 
interaction and it is valid for values of q small enough, 
qR < 1. Eq. JSJ shows that A + B is negative, therefore 
pointing out an instability, in the range sin 2 p < 1/6: 
|0| < cj) c = 24.1° and |tt - <f>\ < (j) c (see note Si). It is 
worth noting that, with increasing R, (f> c keeps constant 
but the maximal value qu ~ 1/R at which (A + B) < 
decreases. The conclusion is that, for any finite value of 
the range R of dipolar interaction, there is a portion of 
the Brillouin zone where oj 2 (q) < 0. This means that the 
ferromagnetic state is unstable for any finite R. 

We have found numerically that the domain wall en- 
ergy Ed w (R), in the case of finite interaction range R, 
has negative values for any R. A limiting value exists 
of Edw(R — * oo) = —0.3570. This result means that do- 
main wall formation is energetically favored for any finite 
R, consistently with the spin-wave instability discussed 
before. 

Note that there is an apparent inconsistency here. 
Above we stated that E^ > 0. How is it possible that 
lima^oo E dw (R) = E dw (w) < and Ef w > 0? The 



4 



answer follows from our previous analysis: the domain 
wall energies E dw (oo) and E^ differ because they corre- 
spond to different ways to take the thermodynamic limit. 
Edw{oo) is evaluated assuming that each spin interacts 
with all spins within a circle of increasing radius R, which 
corresponds to a system of aspect ratio of order one. For 
this system (see the inset to Fig. we find correctly a 
negative domain wall energy. On the other hand, evalu- 
ating E^ with the assumption that each spin interacts 
with an infinite number of other lines (i.e. a system with 
very large aspect ratio), we find a positive domain wall 
energy. A quantitative understanding of the disagree- 
ment between Ed w (oo) and E^ is also possible^ 

The tendency to form domain walls and vortices is very 
sensitive to the presence of local magnetic anisotropies. 
Our final considerations are for the effects of quadratic 
and quartic local anisotropies on spin configurations in a 
dipolc lattice. The anisotropies are defined by 

E am = cos 2 (f>i - K A ^(Sin 4 </>,; + COS 4 (6) 

i i 



where & is the angle formed by spin Si with the x axis 
and K 2 , K4 > so that both single-site anisotropies and 
shape anisotropy favor the x axis. We have evaluated 25 
E anl for a vortex configuration, with the result that 

^(Efr-E&^K^+Kt^ (7) 

which should be added to the vortex dipolar energy (per 
spin), ^(E^-EZ) = %[a(r)-b(r)\nL x ]. For r ~ 
1, the vortex state is favored for L x < c(fl/K), with K = 
max{ii2, K4} and c « 4/10. For large r, a very small 
anisotropy is enough to favor the ferromagnetic state for 
any L x . 

As for future work, it would be interesting to perform 
more sophisticated analyses of the ground state (e.g., 
simulated annealing 24 ), allowing for spin relaxation and 
overcoming of energy barriers. The comprehension of 
temperature effects would be worth studying as well. 
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